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are determined to prove the performance of the proposed method. The 


Image restoration 
Interpolation 
Machine learning 


Noise removal results, when compared with the seven recent state-of-the-art techniques, 
Rectilinear regression show the superiority of the proposed method in terms of visual quality and 
Salt and pepper noise accuracy. The proposed model achieves an average PSNR value of 65.21dB 


and an SSIM value of 0.999 for the reconstruction of medical images, 
proving its accuracy and efficiency. The impulse noise restoration process 
helps the radiologist get a clear visual clarity of the medical image for 
diagnosis purposes. 
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1. INTRODUCTION 

Signal processing is a branch of electronic engineering that employs a variety of processes to 
analyze and modify images, allowing for the most accurate scientific inferences. Digital Image Processing 
algorithms are applied to enhance or reinforce a digital image. The recovery or convalescence of images by 
removing the noise contained in them is one of the most challenging areas of image processing. Pre- 
processing ensures that the unwanted information is isolated from the input data. Image restoration has a vital 
role within the medical imaging field. In medical image processing, the collected medical image data is 
passed through numerous channels or mediums to destinations in several formats, like digital imaging and 
communication in medicine (DICOM). This study eliminates impulse noise from grayscale images using the 
rectilinear regression followed by the linear mean filter. 

Noise recuperation using linear and nonlinear filters has been the subject of extensive research. A 
median filter is a nonlinear filter implemented in an adaptive weighted median filter or a switching median 
filter. On the other hand, the mean is a linear filter implemented for the same by incorporating an iterative 
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mean filter or weighted mean filter. Recent studies show that the combination of median and mean filters is 
applied to the recuperation process. A modified cascade channel [1] is used to eliminate high-thickness salt- 
and-pepper noise is implemented, which makes use of an aggregate trimmer median filter (TMF) and a 
decision-based median filter (DMF). Recent studies exploit machine learning techniques [2], convolutional 
neural network (CNN) models [3] and deep learning techniques [4] using the divide and conquer approach 
for image denoising. In contrast to customary methodologies that learn clean images from uproarious images, 
the model [5] uses residual learning, which takes in commotion from noisy images and afterwards deducts it 
from noisy images in order to get perfect pictures. Garg [6] implements an adaptive weighted minimal and 
maximum median filter, in which the window dimension is increased following the window's corruption-free 
pixel status, and the noisy pixel's value is determined using the weighted median filter. Two novel median 
filter algorithms are introduced in [7], [8] for the recuperation of noise from the grey levels, utilising the rank 
order-based adaptive median filter and the impulse size-based adaptive median filter, respectively. A two- 
stage iterative filter [9] and a non-linear non-local adaptive switching median filter [10] are applied for the 
recuperation process. The non-local mode is most effective when the window has many noisy pixels. A 
hybrid strategy for restoration was applied [11], using the highest frequency non-corrupted pixel in the 
operative region and an applied median filter. A modified pixel median filter for the recuperation process is 
mastered in [12]. The study performed in [13] utilises the frequency of the non-corrupted pixels to compute 
the grey value of the noisy pixels within the region. An adaptive fuzzy detector [14], is applied to detect the 
salt and pepper noise and a weighted mean filter to restore the same. A fuzzy mathematical morphological 
open-close filter [15] uses only uncorrupted pixels for computing the pixel value for eliminating salt and 
pepper noise. A sparse recuperation model for impulse noise is proposed in [16]. 

Image noise restoration has a wide range of real-life applications. Medical imaging is one of the 
relevant areas of image processing in medicine. The medical images generated from various electronic 
devices include noises. For interpretation purposes, preprocessing techniques like visualization, smoothing, 
edge detection, and restoration are applied to the input image. A fuzzy rule-based system for restoration of 
random impulse noise is experimented with in [17]. An S-curve logistic regression technique is mastered in 
[18] for the prediction of diabetes. The linear regression technique [19], [20] is applied for the automatic 
prediction of missing values. A threshold-based technique is employed to reduce noise, and a locally adaptive 
threshold-based edge preservation denoising scheme is used [21]. Digital image transformation [22] 
techniques are applied to remove noise from medical images. A hybrid technique of a modified median filter 
and a Wiener filter [23] was used to restore the noise from images using morphological image processing 
operations. A nonlinear fuzzy filter [24] using a computational intelligence algorithm is applied to medical 
images to remove impulse noise. A type-2 fuzzy kernel for eliminating impulse noise from the grayscale 
image has been proposed in [25], which uses the combination of mean and variance for the restoration 
process. A fuzzy logic technique [26], a direction-weighted mean filter [27], and a Quaternion two-stage 
filter [28] are applied to eliminate noise from medical images. 


2. RESEARCH METHOD 

The block diagram for the proposed scheme is available in Figure 1. The experiment consists of two 
major sections. First, the restoration value for the corrupted pixel is computed by applying the linear 
regression method to each direction. The value obtained in the regression method is given to the mean filter 
in the second stage for determining the recuperating value. If there are no valid directions, then the 
recuperating value is computed using the uncorrupted pixels in the window. A generalised flowchart for the 
proposed model is available in Figure 2. 
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Figure 1. Block diagram 
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Figure 2. Flowchart for the proposed model 


The proposed algorithm for noise restoration has two parts: noise detection and noise restoration. 
Noise restoration algorithm is depticted in Figure 3. The algorithm reads a noisy grayscale or a dicom image 
and applies a border replication algorithm to enhance the image, to perform operations based on a 3x3 
structuring element. For each noisy window detected, it uses the Restore-Window algorithm, explained in 
Algorithm 2 as shown in Figure 4, to compute the restoration value for the corrupted pixel. 


2.1. Hit and miss window 

The restore window algorithm operates on the noisy window to identify all pairs of noise-free pixels 
in all possible directions. If there are no noise-free directions, then a mean filter is applied to the noise-free 
pixels on the window to compute the pixel value for the noisy centre pixel. Otherwise, for each noise-free 
direction, all the pixel values of the noise-free pixels in the window between the two values of the selected 
noise-free direction are used for computing the restoration value using curve fitting linear regression. The 
process is repeated for all noise-free directions, and a mean filter is applied to the computed pixel values to 
get the restoration value. It depends on the nature of the pixels and the origin of the operating region. The 
windows are classified into two categories: miss-window and hit-window. The window in which no 
recuperation process is performed is called the miss-window. The two scenarios leading to miss-window are 
visualized in Figure 5. Figure 5(a) explains the first case in which the origin of the kernel is an uncorrupted 
pixel and no recuperation process is performed. Figure 5(b) depicts the second scenario in which every pixel 
is noisy, including the center pixel, and thus no operation is performed. 
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Algorithm 1 Noise Restoration 


Data: NI = A Noisy grayscale or DICOM Image 
Result: RI = Restored Image 


W + Width( NI) > Width of the Image 
H + Height(NJI) > Height of the Image 
II + Border Replication(N1) > EnhancedInput Image II 
RI + II > Output Image Initialization 


for each, noisypirel € II do 
window + W indow(noisypixel) 
if window is hitwindow then 
pu + Restore — Window( window) > Algorithm 2 
RI; —pv 
end 
end 


Figure 3. Noise restoration algorithm 


Algorithm 2 Restore Window 
Data: IW = A 3x3 window with noisy origin and having at least one corrupted pixels 
Result: pv = Restored pixel value 
DNF + Filter — Noise — Free — Directions(IW) > Filter all pair of pixels in the four directions having 
noise free pixel pair 
NF + All — Noise — Free — Pixels(IW) 
if DNF is ( then 
pu + Mean(NF) > No Direction - Apply Mean Filter 
else 
LReO 
for each, DirectionD € IW do 
min + Minimum — Pizel(D) 
max + Maximum — Pixel(D) 
MLR < Filter — Pixels — Between(min,mar) 
LR + Union(LR,Mean( M LR)) > Direction Based Linear Regression 
end 
pu + Mean(LR) > Restoration Value Computation from the regression values 
end 


Figure 4. Restore window algorithm 


Figure 5. Miss window scenario, (a) noise free center pixel and (b) all pixels are noisy including origin 


2.2. Direction based supervised linear regression 

The window in which the origin pixel is a noisy one with a few noise-free pixels available is called 
the hit window. Based on the number of noise-free pixel directions in the hit window, the recuperation 
process takes place in two approaches. They are direction-based linear regression followed by a mean filter 
and a non-direction-based mean filter. In a 3x3 window, there are four directions concerning the origin pixel, 
where each edge consists of two pixels on the sides, with the noisy pixel in the center. Direction-based linear 
regression is applied only if both the pixels on the edges are noise-free. All the uncorrupted pixel values are 
stored in an array and are given to the linear regression function to predict the best match for substituting the 
corrupted pixel value. The different direction-based linear regression scenarios are classified into tetra, tri, bi, 
and mono-directional. The window consists of only noisy pixels at the origin of the kernel is called tretra 
direction window, shown in Figure 6. All four directions centred on the kernel pixel can be used for 
calculating the recuperation value for the centre pixel as shown in Figure 6(a), and Figure 6(b) shows it in 
shaded color. 
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Figure 6. Tetra direction window, (a) Noisy center pixel with four direction edges and (b) pixel directions in 
different shaded color 


In each of the four directional edges, the linear regression algorithm selects the four most suitable 
values from the list of noise-free pixels. The process is applied only if both the pixels in the directions are 
noise-free. 

D1-{110, LR1, 180} 

D2-{105, LR2, 140} 

D3-{ 100, LR3, 130} 

D4-{200, LR4, 160} 

To calculate the value of LR1, the minimum value is 110 and the maximum value is 180. All non-corrupted 
pixel values between minimum and maximum value, including them are {110, 130,140,160,180}. 

LR1 is the mean of these values. 

LR1=(110 + 130 + 140 + 160 + 180)/S5=144 
Similarly, Other values are calculated as 

LR2=(105 + 110 + 130 +140) /4 = 121 

LR3=(100 + 105 + 110+ 130)/4= 111 

LR4=(160 + 180 + 200) /3 = 180 
Calculate the mean of the values obtained, and that will be the pixel value restored for the noisy pixel. 

PV=(LR1 +LR2+LR3+LR4)/4 

PV=(144 + 121 + 111 + 180)/4 

PV=139 
The calculation procedure explained above is commonly applied to all hit windows having tri-direction, bi- 
direction, and mono-direction. 


2.3. Non direction based linear mean filter operation 
When there is no direction of noise-free pixels available with a pair of noise-free pixels, the 
proposed scheme computes the restoration value by applying a mean filter to all uncorrupted pixels in the 
operating region. Figure 7 demonstrates the scenario. The noise free pixels are {200, 135, 110, and 110}. 
PV = (200 + 135 + 110+ 110)/4 
PV = 139 
The determined value 139 is used to restore the centre pixel. 


Figure 7. Non-direction window 


3. RESULTS AND DISCUSSION 

Experiments with standard images from the literature, sample images from the USC-SIPI dataset, 
and 40 images from the TESTIMAGES dataset are used to assess the proposed system's effectiveness. The 
tool utilized for directing the investigation is MATLAB R2016a. The conventional visual quality estimate 
measures [29] peak signal to noise ratio (PSNR), structural similarity index (SSIM), and Histogram are 
utilised to assess the efficacy of the results produced. The accuracy of the proposed scheme is shown by a 
histogram plot comparison of the original and restored images. The suggested approach restores noise density 
from 10% to 90% in DICOM images with exceptional results. 
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3.1. Peak signal to noise ratio (PSNR) 

PSNR is used to evaluate the reproduction quality of images. The peak signal-to-noise ratio is a 
standard metric for measuring the fidelity between the source and recuperated image. The PSNR ratio is 
determined using (1). 


2552 
PSNR (in dB) = —— py i 
( ) wya Uber Lis (SI (p,a)- RI@.9))? 


(1) 


Where W and H are the image width and height, and SI[p, q], RI[p, q] is the p™ row, q™ column pixel in the 
source and recuperated image, respectively. The prospective framework has experimented with standard test 
images from literature and sample images from the USC-SIPI dataset. The results computed for the PSNR 
values are shown in Table 1. Figure 8 shows the corrupted image and the respective recuperated image for 
the impulse noise density of 30%, 50%, 70%, and 90%, respectively, for the images Pepper, Baboon, and 
USC-SIPI-1 of size 512x512. Figure 9 shows the histogram comparison of Baboon, Family, and USC-SIPI-1 
recuperated images with 90% impulse noise and the source image. 


Table 1. The obtained PSNR values for standard test images and images from the USC-SIPI data set 


Image 10 20 30 40 50 60 70 80 90 
Pepper 44.37 41.02 38.82 37.08 35.54 33.92 31.74 30.02 27.24 
Gold hill 42.32 39.11 37.07 35.43 33.99 32.46 30.68 29.17 26.86 
Girl 44.32 40.91 38.95 37.13 35.33 33.64 31.61 30.16 27.45 
Family 44.07 40.82 38.49 36.61 34.78 33.12 30.76 29.21 26.40 


Baboon 33.60 30.41 2840 26.91 25.66 24.43 23.19 21.93 20.51 
Usc-Sipi-1 42.62 39.31 36.98 35.26 33.40 31.69 30.02 28.29 25.52 
Usc-Sipi-2 41.18 37.87 35.77 34.06 32.53 31.00 29.51 27.90 25.45 
Usc-Sipi-3 42.34 39.19 37.09 35.38 33.83 32.35 30.83 29.27 27.27 


Pepper (70%) sins (90%) 


Pepper (50%) 
ep, ind 


PSNR: 44366 PSNR: 38.820 PSNR: 35.545 PSNR: 31.741 PSNR: 27235 


Baboon (10%) Baboon (30%) Baboon (50%) Baboon (70%) Baboon (90%) 


PSNR: 25.655 PSNR: 20.511 


Usc-Sipi- 1(50%)  Usc-Sipi-1 (70%) nee 1 (90%) 


R STAR 


Ry ee oe 


PSNR: 42.620 PSNR: 36977 PSNR: 33397 PSNR: 30.017 PSNR: 25521 


Figure 8. The noisy Pepper, Baboon, and USC-SIPI-1 images with 10%, 30%, 50%, 70%, and 90% impulse 
noise, and the restored images 


Indonesian J Elec Eng & Comp Sci, Vol. 28, No. 2, November 2022: 777-786 


Indonesian J Elec Eng & Comp Sci 


| 
| 
2s00} 
2000} 
all 
| 
sooo} 


Figure 9. Histogram comparison of the source image with a restored image having 90% impulse noise 
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3.2. Structure similarity index metric 

The structural similarity index metric measures the structural change in the image, which results in 
the degradation of the visual quality of an image with respect to the source image. The SSIM value is 
determined using (2), which takes the source (SI) and recuperated (RI) image as the input. Upon evaluation, 
the computed high SSIM value implies a close similarity between the source and the recuperated image. 
Table 2 shows the SSIM results for the ten standard test images chosen from the literature. 


SSIM(SI, RI) = 


Table 2. The obtained SSIM values for standard test images and images from the USC-SIPI data set 
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sa 


Baboon - 90% Restored 


Family- 90% restored 


Usc-Sipi-1- 90% restored 


Image 10 20 30 40 50 60 70 80 90 
Pepper 0.993 0.985 0.977 0.967 0.955 0.939 0.914 0.880 0.828 
Gold hill 0.989 0.977 0.963 0.947 0.926 0.900 0.859 0.811 0.723 
Girl 0.994 0.987 0.979 0.968 0.954 0.934 0.903 0.869 0.788 
Family 0.995 0.990 0.984 0.974 0.962 0.941 0.908 0.874 0.786 
Baboon 0.977 0.950 0.919 0.883 0.841 0.790 0.726 0.641 0.516 
Usc-Sipi-1 0.989 0.976 0.961 0.944 0.922 0.894 0.860 0.815 0.735 
Usc-Sipi-2 0.986 0.970 0.951 0.929 0.903 0.869 0.826 0.768 0.665 
Usc-Sipi-3 0.983 0.966 0.947 0.924 0.898 0.866 0.826 0.776 0.704 


3.3. Experiment result on TESTIMAGES dataset 

The proposed framework outperformed the majority of recent studies when tested with 40 images 
from the TESTIMAGES [30] data set. Figure 10 shows the comparison of the proposed method with 7 recent 
studies with a noise density of 80%. Figures 10(a) and (b) show the comparison of the PSNR and SSIM value 
of the proposed method with seven recent studies. The results obtained for PSNR and SSIM were compared 
with seven recent studies and found that the proposed framework outperforms all the others with noise 
density from low-density to high-density. Table 3 shows the results of forty images from the TESTIMAGES 
dataset with noise densities of 20%, 40%, 60%, and 80%, and the PSNR and SSIM values obtained after the 


recuperation process. 
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Figure 10. Comparison of the proposed method with 7 recent studies with a noise density of 80% (a) PSNR 
value comparison and (b) SSIM value comparison 


Table 3. Result analysis of 40 images from the TESTIMAGES dataset, with PSNR and SSIM values from the 


proposed approach compared to 7 recent studies 


Algorithm Criterion 20% 40% 60% 80% 
psnr 25.9 22.68 20.57 16.74 

AMF ssim 0.9502 0.9083 0.8401 0.5922 
psnr 29.17 18.28 11.72 RIS 

ACWMF ssim 0.9078 0.5013 0.1493 0.0384 
psnr 32.37 30.37 25.78 20.88 

FDS ssim 0.9587 0.9246 0.8244 0.6401 
psnr 32.15 29.49 27.78 24.72 

BPDF ssim 0.9595 0.9176 0.867 0.7828 
psnr 32.57 29.59 27.2 24.44 

AFMF ssim 0.9611 0.938 0.8938 0.8028 
psnr 34.14 30.3 27.63 24.8 

DAMF ssim 0.9113 0.8752 0.8244 0.7398 
psnr 39.4 35.3 31.701 27.72 

ARmF ssim 0.9844 0.9637 0.9268 0.8498 
psnr 39.2 35.35 31.711 27.99 

Proposed Method ssim 0.9836 0.9638 0.9269 0.8572 


3.4. Image Denoising in Medical DICOM Images 

In medical image processing, one of the standard formats for storing scanned images is called digital 
imaging and communications in medicine (DICOM). The primary reason for using this format is to transmit 
visual data through any communication channel. During transmission, the data may get corrupted by impulse 
noise. This impulse noise in the DICOM image is eliminated using curve fitting and linear regression. 
Figure 11 shows a sample DICOM image with noise restored. The high PSNR and SSIM values prove the 
efficacy of the proposed framework. From Table 4, it is evident that the SSIM value is 100% in most of the 


recuperation processes. 


Table 4. The PSNR and SSIM values for the six dicom images restored from impulse noise 10% to 90% 


Image Criterion 10 20 30 40 50 60 70 80 90 Average 
dicom-1 psnr 77.61 73.96 72.19 70.42 6841 66.32 64.50 61.31 57.62 68.04 
ssim 1.00 1.00 1.00 1.00 1.00 1.00 1.00 1.00 1.00 1.00 

dicom-2 psnr 77.52 74.06 71.98 70.37 68.77 67.28 64.59 62.02 57.14 68.19 
ssim 1.00 1.00 1.00 1.00 1.00 1.00 1.00 1.00 1.00 1.00 

dicom-3 psnr 78.91 75.88 74.03 72.21 70.48 68.84 66.72 64.11 59.78 70.11 
ssim 1.00 1.00 1.00 1.00 1.00 1.00 1.00 1.00 1.00 1.00 

dicom-4 psnr 78.20 74.71 72.97 71.21 69.82 68.03 66.04 63.65 59.92 69.39 
ssim 1.00 1.00 1.00 1.00 1.00 1.00 1.00 1.00 1.00 1.00 

dicom-5 psnr 76.80 73.45 71.40 69.98 68.35 66.53 64.56 62.35 58.88 68.03 
ssim 1.00 1.00 1.00 1.00 1.00 1.00 1.00 1.00 1.00 1.00 
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Figure 11. DICOM image noise restoration from noise densities of 10%, 50%, and 90% 


4. CONCLUSION 

In this study, a combined framework is mastered for the impulse noise recuperation of grayscale and 
DICOM medical images. The developed supervised technique combines the concepts of linear regression and 
linear mean filtering to create a hybrid framework for impulse noise recovery. The experimental result metric 
having high values for PSNR and SSIM proves the standard of the proposed framework. The scheme is also 
applied to remove the noise from medical DICOM images and which helps the radiologist to make inferences 
from the restored image. However, the proposed technique may be further advanced by incorporating deep 
learning techniques into the noise restoration process. 
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